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Abstract 

Multi-soliton solution of the 3-waves problem is represented in ex- 
plicit determinative form. 







1 Introduction 



The problem of three waves in two dimensions arises in different forms in 
many branches of the mathematical physics. For example, this occurs in 
radio physics and nonlinear optics applications, which can be found in [1]. 

Three-wave interaction in plasma, stability criteria, and asymptotic be- 
havior for a general system of three interacting waves, the influence of mu- 
tually different linear damping coefficients on a system of three interacting 
waves of equal signs of energy are described in [2] , where an accurate studying 
waves in plasma were performed. 

Nonlinear Water Wave Interaction and hydrodynamic turbulence includ- 
ing tree wave problem were investigated in [3, 4, 5] 

Furthermore, the three wave problem is also discussed during studying 
nonlinear lattices [7] and photon crystals [8]. 

The present paper must be considered as a direct continuation of the 
previous one [9]. For the convenience of the reader we at first represent nec- 
essary results from [9] marking them into the form of an Appendix, putting 
this on the unequal first place. 



2 Appendix 

2.1 3- wave problem in dimensionless form and its dis- 
crete transformations. 

First off, let us consider the system of equations for 6 unknown functions 
Q,P,A,D.B.E 

P 1 = -QE, A 2 = -BQ, -{Q 1 + Q 2 ) = Q 3 = -PA 

B 1 = -AD, E 2 = -DP, D 3 = -EB (1) 

Using direct calculation it is not difficult to certain that system (1 is 
invariant with respect to the following substitution: T 3 ) 

- _ 1 - B - _ E 
Q ~D' A ~~D' P ~D' 

B = D(^) 2 , E = -D(^) lt ^ = DQ-(\nD) h2 
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The transformation is denoted as T3, which is a discrete transformation. 

Permutating indexes (1, 3) (altogether with corresponding exchanging un- 
known functions) make it possible to obtain 7\ as a discrete transformation 
when the system (1) is invariant. 

S = B(§) 2 , A = -B(^) 3 , | = BP-(lnB) 2 , 3 
And finally the discrete transformation T 2 has the form of: 

A = — , B = Q = , 

D = -£(^) l5 P = £(|) 3 , ^ = EA-(\nE) h3 

In the form presented above, substitutions Tj may be considered as a 
mapping connected to six initial (unbar) functions with six final (bar) ones. 
On the other hand, each substitution should be considered as the infinite 
dimensional chain of equations. For instance, the corresponding chain of 
equations for the case of 71 substitution has the form of: 

on+l r>n jjn An 

7 = -(ln£ n ) 2 3, D n+1 = B n ( — ) 2 , A n+1 = -B n ( — ) 3 (2) 

ptn An 
B n ' ** £n 

In the first row, we have the lattice as if it were system connected to 3 
unknown functions (B, D, A) in each point of the lattice. The first chain for 
B functions is definitely well known as a two dimensional Toda lattice. 

2.2 Some properties of the discrete transformations. 

All of the discrete transformations constructed above are invertable. This 
means that an unbar unknown function may be presented in terms of the 
bar ones. For instance, T 3 _1 looks as: 

D = , B = — E = , 

Q Q Q 
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P=-Q(Z)2, A = Q(~) U 2 = 5Q-(ln<5)i >2 

Consequently, it is not difficult to certain by direct computation that 
discrete transformations Tj are mutual commutative {T{Tj = TfTA in the 
solutions of the system (1). Moreover, others, such as T{F 2 = T 2 T\ = t 3 , can 
be also handled in the same way. 

Thus, from each given initial solution Wq = {A, P, Q, E, B, D) of the 
system (1), it is possible to obtain the chain of solutions labeled by two 
natural numbers (li, l 2 , or (7 3 )) and the number of applications of the discrete 
transformations (Ti, T 2 , T 3 ) associated to it (as it was noticed above, T{T 2 = 
T 2 T 1 = T 3 ). 

Arising chain of equations with respect to (D, B, E) functions are de- 
finetely two-dimensional Toda lattices. Their general solutions in the case 
of two fixed ends are well-known [13]. As the reader will see, soon this fact 
allows constructing the many soliton solutions of the 3-wave problem into 
the most straightforward ways. 

2.3 Resolving of discrete transformation chains. 
2.3.1 Two identities of Yacobi 

The first Jacobi identity: 

(T n -ia l \ /T„_ia 2 \ /T„_ia 2 \ (T n -ia l \ 

Un \ &Vn J Un I b 2 r 22 J ° n \ Vt 12 ) ° n I b 2 r 21 ) ~ 

£> n _i(r n _i)D n+ i b l T 1X T 12 

V b 2 T 2 iT 22 J 

where a\b l are (n — 2) dimensional columns (rows) vectors, r^j components 
of 2-th dimensional matrix. 
The second Jacobi identity: 

Dn ( ) D n+1 I fop I - D n J D n+1 I foiM I = 
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D r , 



T n -\a l 



( T n _ia l a 2 



D 



n+l 



b 2 pr 
\ b l Ta 

These identities can be generalized in the case of arbitrary semi-simple 
group. The reader can find these results in [11]. 



2.3.2 Resolving of the discrete lattices 

Let us take an initial solution in the form of: 

Q = A = P = 0, B = B(2), E = E(l). D 3 = -BE 



(3) 



Application to this solution to each of the inverse transformations Tf 1 is 
meaningless because of arising zeroes in denominators. The chain of equa- 
tions under such boundary conditions is what we call the chain with the fixed 
end from the left (from one side). 

The result of applications to such initial solution by 1$ times T3 transfor- 
mation looks as (in order for checking this fact only two Yacobi identities of 
the previous subsection are necessary). 



-iy 



-i^3zi 5 D (h = (_i)'»^2±I, A = l 



*3 



A B 

A(h = (-lya-ia, 

A «3 



A B 
^0 



A i3 ' 



A E 
A / 3 



^0 



A_i = 0. 



(4) 



where A„ are minors of the n-th order of infinite dimensional matrix 



( D D 2 D 22 
Di D12 D 12 2 
D\\ T)\\i -D1122 

v 



(5) 



/ 



and A B ,A B are the minors of l 3 order in the matrices of which the last 
column (or row) is exchanged on the derivatives of the corresponding order 
on argument 1 of E function (on argument 2 of the B function in the second 
case) . 
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In the following next notations will be used. W l:3 ' 11 , (W l:i,l2 )~ this is the 
result of application of discrete transformation T l:i T l1 (T l:i T l2 ) to the cor- 
responding component of the 3-wave field. A^ ^ (A; 3 / 2 ) - determinant of 
^3 + h (h + h) orders, with the following structure of its determinant matrix. 
The first I3 rows (columns) of it coinside with matrix of (5) and last l±, (I2) 
rows (columns)are consructed from the derivatives of B, (E) functions with 
respect to arguments 2, (1). 

The result of additional application to l\ times Ti transformation to the 
solution (4) looks as : 

ph,h — ^h,h-i r>l 3 ,h _ \w a - 1 A ; 3 -1 — A E 

Qh,h _ ( iy-i+h-1 ^h-l,h jjL z ,h _ ^ ^y-s+h ^la+hh ^ 

A Wi ' Ai 3 ,h 

£jh,h = ( iy 3 +h ^ h+l,h-l ^ ^l 3 ,h = ( iy 3 +h ^3-l,h+l ^ 

A /3,/l A «3,h 

We do not present the explicit form for components which can be 

obtained without any difficulties from (6) by corresponding exchaging of the 
arguments and unknown functions. 



3 Multi-soliton solution of the 3-waves prob- 
lem 

3.1 General consideration 

The system (1) allows reducing (under additional asumption that all opera- 
tors of differentiation are the real ones d a = d* a ) . 

P = B\ A = E\ Q = D* (7) 

In this case the system (1) is reduced to three equations. 

B 1 = -DE*, E 2 = -DB*, D 3 = -BE (8) 

for three complex unknown functions (E, B, D). This system of equations is 
definetely 3-wave problem. 
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In [9] it was shown that begining with an initial solution (3) after the 
(h, hih)) steps of discrete transformation leads to a solution of the reduced 
system (8) if the following conditions are satisfied: 

^2Z 3 +l,2h = ^2*3+1,2*1-1 = A 2 * 3) 2* 1 +l — 

^2*3-1,2*1+1 _ ^_-Q* 3 +*i£>* ^2*3-1,2*1 _ (_iyi-l£)* ^2*3,2*!-! _ (_-Q*3£* 
^2*3,2*i ^2*3,2*i ^2*3,2*! 

(9) 

Resolving the first line of above equations is as follows: 

2*3 

D = ]T/W(2), 
k=l 

To resolve the second line ones the following parametrization is sufficient: 

2*3 2*3+2*1 

s=l 5=1 
2*3 2*3+2*1 

£ = ]Te s e A *\ B= b s^ s2 

s=l 5=1 

All numerical parameters c^,dg,e s ,bs are connected by equations (9) and 
one additional equation connects with initial conditions (3). 
The initial condition leads to: 

2*3 2*3 ( r -l\k 

fc=l «=i A « + 

and finally we obtain for D : 

2*3,2(* 3 +*i) A,i ft „ s2 
s =l,5=l A * + MS 

Thus, solution to be determined are defined by the pair of (2/ 3 ) parmeters 
e s ,X s and 2(l 3 + li) pairs parameters bs,^s- Equations (9) give some ad- 
ditional limitations (of reality) for these parameters. Below, we find these 
limitations and then I3, l\ is given by this explicit formula (6), in which con- 
ditions of reducing are satisfied. 



6 



3.2 Determinant camputation. The case of denumer- 
ator. 

At first, let us calculate the common determinant (see (9)) -A^^ using 
explicit expressions for initial functions D, E, B. The determinant matrix of 
dimension 2(l 3 + l 1 ) x 2(7 3 + Zx) may be presented as product of two matrices. 
The first one has the block form 2/ 3 x 2/ 3 matrix L(e, A; 1) in upper left 
angle and zeroes on other places except of the unity on main diagonal. The 
elements of the matrix L are the following = A^e^e^ 1 . The result is 
obvious: 

AS, >2Il = Det(L) = ri^ 21s (A 1 ,...A 2i3 )e(2 1 W (11) 

3 = 1 

where If is a Vandermond determinant. The second matrix in its turn can 
be represented as product of another two ones. The first 2l 3 lines of the 
first matrix have the matrix elements t-4 — (s number of line, S number of 
column). The remaining 2/i lines of this matrix are the usual Vandermond 
matrices W^g = l^s' 1 - Computation of this determinant is not complicate 
problem with the following result: 

s=2l 3 ,S=2(l 3 +h) j 

A 2i S) 2ii = W 2h{^ ■■Mh)W 2i i z+h) {n 1 , ...^ 2{h+h) ) [] — — 

3=1,5=1 A « + 

Finally, the matrix elements of last 2(7 3 + li) x 2(/ 3 + li) matrix have the 
following analitical structure ^e^/^- 1 And in consequence: 

A l,2 h = ^ bsW 2(h+h) (^ ...^ (i3+!l ))e ( SS +,l) "^ 2 

5=1 

Thus, for A 2 i 3i 2/i summating all results above we obtain: 

s=2l 3 ,S=2(l 3 +h) ^ 2l 3 2(l 3 +h) 

^ 3 (A 1 ,...A 2 ,3)w 2 2 { , 3+ , l) ( / x 1 ,... / / 2( ^ 1) ) n y— — n^ eAsl n ^& 

s=l,5=l As ~+~ s=l 5=1 

(12) 

3.3 Determinant computation. The case of B function 

In the process of computation of A 2 ; 3j2 / 1 _ 1 , the determinant (2(/ 3 + li) — 1) x 
(2(Z 3 + li) — 1) matrix may be represented as product of two matrices. The 
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first one coincides with matrix L of the previous subsection. Determinant of 
this matrix has been already calculated in (11). Second matrix in its turn 
may be represented in the form of the product of two rectangular matrices 
of dimentions (2(Z 3 + h) - 1) x 2(Z 3 + Zi)) and 2(Z 3 + h) x (2(Z 3 + Z x ) - 1), 
correspondently. The structure of these matrices are described in the previous 
subsection. The determinant of the degree 2(Z 3 + l±) — 1 is equal to sum 
of products of all minors of 2(Z 3 + Zi) — 1 orders of both matrices. These 
determinants were computed also in the previous subsection. The finally 
result is as follows: 

2(h+h 2(« 3 +/i) s=2l 3 ,S=2(l 3 +h) j 

w 2i 3 (^i, •••A2/ 3 )M / 2 2 (/3+ii) _ 1 (/xi, ..^ 2 (h+h)) n ^e^ 2 n — 

i=l 5=1, S+i s=l,S=l,& As + ^ 5 

Taking into account the last equality from (9) and substituting in it all the 
results for determinant calculations above, we obtain: 

2(l 3 +h) 2l 3 2(l 3 +h) 1 2(l 3 +h) 

E U(*. + Hs) II (^-^)- 2 e-^ 2 - = (-l)^ E bJe^S 2 

5=1 s=l K=\,K=£S ° s 5=1 

The last equality can be satisfied only under condition: 

fi*s = -^ PS , S^(PS) 

where P operator of permutation of 2(Z 3 +Zi) numbers with obvious propertity 
P 2 = 1. Comparing terms under the same exponents, we obtain: 

2l 3 2(l 3 +h) 

(-l) l3 b*ps = Il( X s + Vs) II {»K-»s)- 2 e-» s2 - 

s=l K=1,K^S °S 

This is typical for soliton theory connection between the "energy" and am- 
plitudes. 

3.4 Determinant computation. The case of E function. 

Numerical parameters c s are connected to the first equation (9). To find 
these relations it is necessary to calculate A 2 / 3 -i i 2« 1 +i- The technique of 
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these calculations are the same as in previous subsection and finally equality 
(9) for E functions looks as: 

2l 3 2l z 1 2{h+h) 2/ 3 

E^^E-e-^ 1 n (*s + »s) IT (A,-A s )- 2 

s=l s=l e s S=l k=l,k^s 

As in previous subsection, the equality above restricted values of parameters 
A as: 

K = — Ap s 

where P is some permutation of the group of 2l 3 numbers satisfying the 
condition P 2 = 1. Comparison the coefficients at the same exponents in both 
sides of the last equality leads to the following conditions on parameters e s : 

2(l 3 +h) 2l 3 1 

<= n (Xps+^s) n (A fe -A s )^ 2 — 

5=1 k=l,k^s C ( Ps ) 

3.5 Determinant computation. The case of D function. 

Calculation of the determinant defining D function and corresponding condi- 
tion of reality (9) does not add new limitations to the numerical parameters 
obtained above in the last two subsections. 

4 Outlook 

No one of the authors is a specialist in the field of applications of 3-wave 
problem to certain physical phenomenon. Due to this reason, we can not 
discuss the usefulness of such applications. 

From a computational point of view, it seems that here we have the first 
paper having a systematical investigation of a multicomponent integrable 
system. It is interesting to compare the results in the case of multicom- 
ponent (the simplest one) integrable system connected with the algebra A 2 
and numerous integrable systems connected with A\ algebra. The strategy 
of computations in both cases is basically the same. Beginning from some 
simple solution connected with the upper triangular nilpotent subalgebra 
((P = A = Q = 0) in the case of the present paper after corresponding 
2n steps of discrete transformation, we come to solution connected with the 
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lower triangular nilpotent subalgebra. Connection with hermitianity of these 
solutions leads to limitation of (arbitrary up to now) numerical parameters 
of the problem. In the case of two component integrable system, for instance 
in nonlinear Scrodinger equation, this dependence looks as [?]. 

2i 3 1 

<= n (^-A s r 2 — , \: = -\ Ps 

k=l,k^s C ( Ps ) 

In the case under consideration, this dependence modified by the formulas 
of subsections B, E. This comparison offers possibility assuming that each 
simple root of semisimple algebra borns its own systems of pairs of parameters 
amplitude-phase and numbers of these pairs are arbitrary and only by them 
multi-soliton solutions are defined in the case of multicomponent integrable 
systems. 

In our consideration, it was unable to connect constructed solutions with 
L-A pair formalism and authors have no idea how this can be done (if pos- 
sible). 
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